Controllability and universal three-qubit quantum computation with trapped electron 

states 
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We show how to control and perform universal three-qubit quantum computation with trapped 
electron quantum states. The three qubits are the electron spin, and the first two quantum states of 
the cyclotron and axial harmonic oscillators. We explicitly show how universal 3-qubit gates can be 
performed. As an example of a quantum algorithm, we outline the implementation of the 3-qubit 
Deutsch-Jozsa algorithm in this system. 
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I. INTRODUCTION 



Since the eighties, much effort has been put into the 
study of quantum computers, and various proposals for 
the physical implementation have been put forward in 
very different fields such as those of cold trapped ions 
[H, quantum dots 0, NMR 3 and neutral atoms 
Experimental implementations have, however, been few 
and limited to a low number of qubits. It is therefore 
interesting to study alternative candidates for quantum 
computers. 

In the present work we study a scheme put forward 
for instance in ^5^ in which an electron trapped in a Pen- 
ning trap is considered a possible realization of quantum 
logic. Advantages such as an almost complete absence 
of decoherence mechanisms and already obtained good 
experimental accuracy suggests this to be an interesting 
direction of study. We show that it is possible to perform 
universal three-qubit quantum computation in this sys- 
tem, and outline the implementation of the three-qubit 
Deutsch-Jozsa algorithm in this system. 

In our scheme, quantum information is stored in an in- 
ternal degree of freedom (spin) and two external degrees 
of freedom (cyclotron and axial motion) of the trapped 
electron. For two qubits, proposals have already been 
put forward in 0] and p. Here, instead of restricting 
ourselves to two qubits, we exploit an additional degree 
of freedom in a single electron and consider three qubits. 

Note that although universal quantum computing for 
three qubits of a single electron has already been demon- 
strated theoretically in 7:], this proposal relied on small 
relativistic effects which lead to anharmonicities for the 
cyclotron oscillator. In our scheme, we treat the cy- 
clotron oscillator as essentially harmonic (to experimen- 
tal resolution) and we use a traveling field that is exper- 
imentally well-established Q. This scheme is distinctly 
different from schemes in alternative trap configurations 



such as in Refs. @ and where a linear array of Pen- 
ning traps was investigated, and in [lH, where an array 
of planar Penning traps was considered. 

In the following section, we briefiy describe the phys- 
ical system and the theoretical model that describes it. 
Then we describe the control mechanisms and show the 
control equations. Section 4 demonstrates that univer- 
sal 3-qubit quantum computation is possible using these 
controls. Subsequently, we present an explicit example of 
the 3-qubit Deutsch-Jozsa algorithm that can be imple- 
mented. In the last section, we discuss the experimental 
challenges to the implementation of this scheme. 



II. THE SYSTEM 

The physical system consisting of a single electron 
trapped in a Penning trap has already been described in 
some detail in for instance p^. To summarize, the elec- 
tron experiences a magnetic field B = Bz and a static 
quadrupole potential 



V^Vo 



2<P 



where Vb is the potential between the trap electrodes 
and d is a characteristic length of the trap p^]. The 
Hamiltonian, given by 



ifo = 7^ (P - eAo)' 
Zra 



eV ^JL■B 



where Aq = x r and jjb — j^u, can be recast in terms 
of three harmonic oscillators [131 
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namely the cyclotron, axial and magnetron motions. 
The frequencies are given by = [uc + Cjc)/'!, ujm — 

(CO, - ^,)/2, uz = y^, = V^^,-2ul = ^ 

and LUs — ^^2m ■ ^P^'^ ^^'^ axial degrees of 
freedom are used as qubits whereas in spin and the 
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cyclotron oscillator constitute the qubits. Here we use 
the spin of the electron and the axial and the cyclotron 
oscillator states for the three qubits. The logical states 
|0) and |1) corresponds to || ) and jt ) for the spin and the 
Fock states |0) and |1) for the cyclotron and the axial os- 
cillators. The states are denoted by \jnl) where j, n, I are 
the logical states for the spin, cyclotron and axial states, 
respectively. The frequencies corresponding to the qubit- 
energy splittings are of the order GHz, GHz and MHz, 
respectively. See Fig. 1 for an energy-level diagram. 



III. CONTROLLING THE SYSTEM 

The spin qubit is the only qubit which can be controlled 
directly. This is done using a small transverse magnetic 
field 12] 

h{t) — b (cos((jjt + (j))x + sm{u!t + (t))y) 

If uj is close to LUs this leads to the following Hamiltonian 
in the interaction picture 



9|e|fc 



- (J- exp(i0)) 



where Q. = In the basis | Inl), \ ]nl) this interaction 

implies the evolution matrix 



cos(f ) 



-ie'"^ sin(^) 
«'^srn(f) cos(f) 



subsequently referred to as a,ps{d, (j)) pulse, where 9 — fit. 

The cyclotron and axial qubits cannot be directly con- 
trolled. Cyclotron transitions can, admittedly, be ad- 
dressed by setting up a proper vector potential but 
since the levels are equally spaced, population will in- 
evitably leave the computational subspace [31. To con- 
trol the spin-axial transition a traveling magnetic field is 
set up 

h{t) ~ b (cos(wi + (f> — kz)x -f sm{ix}t + 4) — kz)y) 



which leads to the following Hamiltonian in the interac- 
tion picture 

H'sa =^{<y+ exp(»(r7(4e-^* + a,e— - A • i - </>)) 
+ exp(-i(r;(4e^"'^* -t- a^e^^'^^*) - A • i - 0))) 



where £7 — 4^^, A 



„ , ^ — cj — Ws and rt = k\ tt^ — . For 

Lu — ujsa — — only the levels | | n[l + 1)) and 
I "[nl) are connected. In the basis | [n{l + 1)), | Inl) the 
evolution matrix is given by 



cos(|) 



sin(|)^ 



l^e-*^sin(|) cos(|) 



(1) 



where 0; = fit • 776 "^^ /'^ ^ L\{rf). (a;) is a general- 
ized Laguerre polynomial. In the following this evolution 
matrix will be referred to as a Psai^^ 4') pulse. Note that 
this does not assume the Lamb-Dicke limit, so the extent 
of the zero-point motion of the harmonic oscillator can 
be a significant fraction of the wavelength of the field. 

A spin-cyclotron interaction can be set up using a mag- 
netic field as described in |6j 

h{t) = b{xx A- yy) cos{ijjt + (j)) 

This leads to the Hamiltonian 



E uj — ujsc — ^s~^'c only the levels \inl), W {n — 1)1) are 
connected. In the basis |in/), |t (^t- — 1)0 the following 
evolution operator is thus obtained 



^ / cos(^) le^^sH^) 
Ue-^4> sin( V^) cos(^) ^ 



(2) 
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FIG. 1: Energy-level diagram. States are denoted |jnZ) , where 
j is the spin state, n the cyclotron state and I the axial state. 
Levels marked with bold are included in the computational 
subspace. Solid, dashed and dotted lines indicate spin, spin- 
axial and spin-cyclotron transitions, respectively. 



where 9 — ~9I^Bby ^mhu ^' the following this evolu- 
tion matrix will be referred to as a Psc{9, 0) pulse. 

Fig. 1 shows that the three fields connect all the eigen- 
states of the spin-axial-cyclotron system. In addition, 
the levels are connected in such a way that the system is 
eigenstate controllable [31 . That is, population can be 
coherently transferred from any eigenstate to any other 
eigenstate. For example, consider the set of eigenstates 
illustrated in Fig. 1. The condition for eigenstate con- 
trollability is that the pulses of frequency LOsa, and oj^c 
must be applied sequentially^ and not simultaneously. For 
example, let us say we want to transfer the lOOO) state to 
the 1111) state. We can do so by the sequence: Ps(7r, 0i), 
Psa(7r, 02), Ps(7r, 03), psc(7r, 04), Ps(7r, 05), where 0i, 02, 
03, 04 and 05 are arbitrary phases. This is similar in 
spirit to the eigenstate transfer schemes for trapped ions 

m. 



3 



IV. UNIVERSAL QUANTUM COMPUTATION 

A universal set of gates consists of the Hadamard, the 
T-gate and the controUed-NOT gates [l3|: 



p,,(7r,0), p,c(V\/2,^/2), Psc{n,0), (VA 7r/2). To 
convert it to a standard CZ gate 




CNOT 



In the following, we will demonstrate how to implement a 
universal set of gates for the current system by applying 
the interactions described above. 

For the spin qubit, the single-qubit gates are easily im- 
plemented. Thus the T gate is performed up to a global 
phase factor using the two pulses Ps(7r, 7r/8), Ps(7r, 0). 
An arbitrary phase gate ( g gi**" ) is in fact performed just 
by Ps(7r, (p/2), ps(7r, 0). The Hadamard gate is obtained 
from the pulse sequence Ps(7r, — tt), ps(7r/2, 7r/2). 

Implementing a swap gate 



SWAP 



/I 0^ 

10 

10 

Vo 1/ 



between the spin and cyclotron qubits is complicated by 
the fact that driving the transitions |01^) «-> |10Z) simul- 
taneously drives the transitions |02^) -f-> |11^) and popula- 
tion thus leaves the computational subspace. This prob- 
lem can, however, be overcome using composite pulses 
p^ . a method that has already been applied to for in- 
stance trapped ions [l^ and trapped electrons [g']. Ap- 
plying the pulse sequence Psc(7r/%/2, 0), Psc(27r/\/2, 0^), 
Psc(7''/\/2, 0), where (j)s ~ arccos(cot^(7r/\/2)), imple- 
ments the swapping gate up to a phase factor of e*" 
for |01) and -e"'" for |10), where a « -0.8652. To 
convert it to a proper swap gate we suggest the follow- 
ing. First, ps(7'', a/2), Ps(7r,0) followed by the swap- 
ping gate and subsequently the pulses Ps(t^, (tt — a)/ 2), 
Ps(7r,0) implements a swap gate up to a phase of tt 
for 1 00). As demonstrated in @ the pulse sequence 
Psc(7r, 0), Psc{'n/V2, 7r/2), Psc(7r, 0), Psc(7r/V2, it/2) leads 
to a controUed-phase gate endowing all states but |00) 
with a phase of tt. Thus the phase of |00) is corrected by 
applying this controUed-phase gate. 

Now implementing any singie-qubit gate for the cy- 
clotron qubit is easy, since IxS ^ SWAP -{Sx I)- SWAP , 
where S denotes an arbitrary single-qubit gate. Imple- 
menting for instance a T gate, however, simplifies since 
sandwiching a T gate for the spin qubit between the 
swapping sequence and the same sequence with all phases 
offset by tt leads to a T gate for the cyclotron qubit. 

To implement a controUed-NOT gate between the 
spin and the cyclotron qubit, we first implement the 
controUed-phase gate as before by the pulse sequence 



a 

10 

10 

vO -ly 



we simply sandwich it between NOT gates on both 
qubits. Finally, applying a Hadamard gate to the 
cyclotron qubit before and after the CZ gate leads to a 
CNOT gate. 

To implement a swap gate between the spin and the 
axial qubits Psa pulses are naturally employed. As before, 
it is a problem that driving for instance the transitions 
|Onl) <-> |lnO) simultaneously drives the transitions 
|0n2) <-> |lril) and population thus leaves the compu- 
tational subspace. Noting that 9i=i = ^(2 — ff)9i^Q 
we can, however, in the limit where ?/ w (for which 
di=i ~ \/2di=o) use the composite pulse idea employed 
for the spin-cyclotron interaction. If all phases are offset 
by —tt/2 (cf. formulas ([T]) and ([2])) we can use exactly 
the same sequences as for the spin-cyclotron case to 
implement the desired gates. Notice that this method 
can also be applied for a nonzero rj. For instance, 
choosing rj = 2 leads to — — \/2^;=o and the pulses 
above are immediately applicable. 

A controUed-NOT gate between the cyclotron and 
the axial qubits is implemented by the sequence 

C NOT^yQio spiri} C N OT spiji fixial 1 C N OT^yQiQ ^p-iyi^ 
C NOTgpj^ji fixial ■ 

So we have now demonstrated how to perform univer- 
sal computation with three-qubits in a trapped electron 
scheme. 



V. IMPLEMENTING THE DEUTSCH-JOZSA 
ALGORITHM 

In this section, we will demonstrate how to implement 
the 3-qubit Deutsch-Jozsa (DJ) algorithm [13]. The re- 
fined version suggested in [2l| allows us to test functions 
for three qubits instead of just two with three physi- 
cal qubits available, as the control register is eliminated. 
Testing functions on three qubits compared to two is in- 
teresting, since three is the minimum number of qubits 
required to solve the D J algorithm in a non-classical way 
[2]| . The algorithm works in the following way: 1) Ini- 
tialize the system to |000), 2) Apply a Hadamard gate 
to each qubit, 3) Apply the unitary operator Uf defined 
by Uf\x) — (— l)''''^)|a;), 4) Apply a Hadamard gate to 
each qubit, 5) Measure the state of the qubits. If the 
measurement returns |000) the function is constant, else 
it is balanced. 

We have already demonstrated how to implement a 
Hadamard gate for each qubit, so we only need to im- 
plement the Uf gates. Now, for three qubits there are 
2 constant functions and 70 balanced functions. These 
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functions are indexed as in [24I using their outputs ex- 
pressed as hexadecimal numbers such that for instance 
the function with the output /(0)---/(7) = 00001111 
is denoted /of- Since Uf^ — —Ufpp__^ we only need to 
implement 1 constant and 35 balanced functions. The 
constant function corresponding to /oo is easily imple- 
mented since Uf^^ is equal to the identity. To implement 
the balanced functions we refer to [12], where NMR se- 
quences for the 35 balanced functions are written down. 
From this work it follows that all functions can be con- 
structed from gates on the form 



1 

-1 



1 

±i 



/I ON 

±i 

±i 

VO 1/ 



The single-qubit gates are easily implementable for the 
spin qubit. The first one, corresponding to /^(Tr) in 
[22| . is obtained from the sequence Ps{n, Tr/2), ps(7r, 0), 
whereas the second one, corresponding to /2(±7r/2), is 
given by Ps(7r, ±7r/4), Ps{tt, 0). Using SWAP gates these 
gates are transformed to single-qubit gates for the cy- 
clotron and the axial qubits and thus all necessary single- 
qubit gates listed in are obtained. Likewise the two- 
qubit operations are easily achieved. Thus the sequence 
CNOT,{l ) , CNOT implements the desired two-qubit 
gates, that is, the gates Jij{±n/2) in [22*1. It follows that 
all the Uf functions are implementable, and therefore 
that the Deutsch-Jozsa algorithm can be executed in our 
system. 



VI. DISCUSSION AND CONCLUSION 

In conclusion, we have demonstrated universal quantum 
computation for trapped electron three-qubit states. We 
have explicitly determined a universal set of gates by find- 
ing the necessary pulse sequences. We have also demon- 
strated that the three-qubit Deutsch-Jozsa algorithm can 
be implemented in this system. 

Thus, this work expands current work in the area of 
trapped-electron quantum computing to exploit another 
degree of freedom in each electron. Besides constituting 
an interesting system in itself in that it can be used to 



test for example the Deutsch-Jozsa algorithm, this work 
could also be useful for studying many-particle gates via 
the Coulomb interaction between electrons in planar Pen- 
ning traps (the system is not scalable with the trap con- 
sidered in this proposal). Each particle could have three 
qubit degrees of freedom and one of these could be used 
for communication. Such a situation is highly interesting 
for the case of error correction. 

However, there are significant experimental challenges 
to the implementation of this proposal at present. The 
primary challenge is that the axial motion cannot yet be 
cooled to the ground state (which is necessary for initial- 
ization). There are proposals to do so using the cyclotron 
motion, but these are still in their early stages. Neverthe- 
less we expect that our work will motivate rapid experi- 
mental progress in this direction, much as early quantum 
computing proposals [ij motivated the motional cooling 
of trapped ions [23l| about a decade ago. A second chal- 
lenge is that all the three-qubit states cannot be mea- 
sured directly. The spin qubit can be measured dl , as 
can the axial qubit with a low temperature circuit [24| . 
For the cyclotron qubit, a scheme adopted in [l3] can 
be used. In that work, in which the spin and cyclotron 
levels are qubits, it is demonstrated that the spin and cy- 
clotron levels can be coupled with the axial motion and 
subsequently the axial motion can be measured to get 
a cyclotron qubit measurement. Although this system 
has long coherence times, there are notable sources of 
decoherence that set a limit on the number of gate op- 
erations that can be performed, for instance fluctuating 
magnetic and electric fields, thermal noise and intensity 
noise from the microwave sources. In spite of these exper- 
imental challenges to be overcome, this system provides 
an attractive quantum computing paradigm to be inves- 
tigated. 
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